Enhanced Microlensing by Stars Around the Black Hole in the Galactic Center 

Tal Alexander 
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218 

and 

Abraham Loeb 

<^ ■ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

O' 

The effect of stars on the lensing properties of the supermassive black hole in the Galactic 

Center is similar to the effect of planets on microlensing by a star. We show that the dense stellar 

cluster around SgrA* increases by factors of a few the probability of high-magnification lensing 

events of a distant background source by the black hole. Conversely, the gravitational shear of the 

black hole changes and enhances the microlensing properties of the individual stars. The effect is 

largest when the source image lies near the Einstein radius of the black hole (1.75" ± 0.20" for a 



CO 

> 

Tj" \ source at infinity) . We estimate that the probability of observing at least one distant background 

^-r ' star which is magnified by a factor > 5 in any infrared snapshot of the inner ~2" of the Galactic 

OS ■ Center is ~ 1% with a .ff-band detection threshold of 20 mag. The largest source of uncertainty 

in this estimate is the luminosity function of the background stars. The gravitational shear of the 
black hole lengthens the duration of high-magnification events near the Einstein radius up to a 
few months, and introduces a large variety of lightcurve shapes that are different from those of 
isolated microlenses. Identification of such events by image subtraction can be used to probe the 
mass function, density and velocity distributions of faint stars near the black hole, which are not 
"^ ' detectable otherwise. 
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1. Introduction 

Deep infrared observations of the innermost region around the massive black hole (BH) in 



the Galactic Center (GC) reveal numerous point sources (Genzel et al. 1997; Ghez et al. |199S ). 
Most of these are probably stars orbiting deep in the potential of the BH. However, the BH will 
also gravitationally lens and magnify any background star that happens to lie behind it, and so 
a small fraction of these point sources could be lensed images of distant background stars. The 
possible existence of such images in the innermost GC, beyond being a probe of the BH potential, 
may also affect the apparent star counts, radial stellar density distribution and infrared luminosity 



function. Previous investigations of gravitational lensing by the BH (Wardle & Yusef-Zadeh 1992 



Alexander & Sternberg 1999| ) suggested that lensing effects do not play a major role in present-day 



observations of the GC. However, the BH is surrounded by a very dense stellar cluster (e.g. Genzel 



et al. 1997] ), whose contribution to the lensing by the BH has so far been neglected. Because the 
stellar mass is not smoothly distributed around the BH but is composed of discrete point masses, 
its effect on the lensing properties of the BH is much larger than one may naively estimate by 
adding the stellar mass to that of the BH. 

The effect of a star on lensing by the BH is analogous to the effect of a planet on microlensing 
by its star. The latter problem has been studied in great detail in the context of microlensing 



searches for planets (Mao & Paczynski |l99l| ; Gould & Loeb |1992| ; Bolatto & Falco |1994| Gaudi 
& Gould |l99ll Wambsganss 1997; Peale |i"997| ; Griest & Safizadeh [l998| ; Gaudi, Naber, & 
Sackett 199J; Gaudi & Sackett 1997 ). In particular, Gould &: Loeb ( 1992| ) first showed that the 
cross-section for magnification by the planet can be increased by up to an order of magnitude due 
to the gravitational shear of the star; the effect being most pronounced when the planet lies near 
the Einstein radius of the star. Aside from the change in scales, the same result should apply to 
SgrA*. In the GC, the primary massive lens is the BH and the secondary low-mass lenses are the 
stars around it. If a star happens to pass near the image of a source that is being lensed by the 
BH, then the magnification of the image could be enhanced significantly. The BH shear increases 
the probability for high magnification events by individual stars. Miralda-Escude & Gould (2000) 
pointed out that the existence of a cluster of stellar mass BHs around the massive BH could be 
detected by this effect. In this paper we consider normal stellar lenses and calculate the probability 
of high-magnification events in the stars-BH system in comparison to the naked BH case. 

The statistics of microlensing by stars in the GC region is very different from that in 
the Galactic disk or halo. In the latter case, the optical depth for microlensing is very small, 
r ~ 10~ 7 -10 -5 , but there are many background sources (see, e.g. review by Paczynski 1996). In 



contrast, the optical depth for microlensing by stars within the inner arcsecond of the GC is as 
high as ~0.1 due to the high stellar density there (see §|j), but it is not clear whether there is 
a sufficient number of bright background sources. In fact, due to their large distances and high 
extinction by dust, most sources behind the GC will only be observable while being magnified 
during a lensing event. The microlensing events are transient; any distant source behind the GC 
will be microlensed for a fraction of the time due to foreground stars which are passing in front of 
its images. 

This paper is organized as follows. The method of our calculation is described in §g; the 
models of the stellar lenses and background stellar sources are described in §|3|; and the numerical 
results are presented in §Q. We discuss the main implications of our results in qq. 



Method 



We adopt the semi-analytic method of Gould & Loeb ( 1992 ; see also the original discussion 
by Chang &; Refsdal 1979] , 1984 ) to describe the effect of a low-mass secondary lens (a star) on 
the image of a point source produced by a high-mass primary lens (the BH). The limitations of 



the point source assumption are discussed in §||. We consider rare, high-magnification events for 
which the possibility that more than one star affects the microlensing event at any given time can 
be neglected. The original, unperturbed image position at x im is displaced by ^ = (£«,??«) due to 
the perturbing star at a position £ p = (£ p ,%>) relative to the unperturbed image (Fig. [[]). We 
express the source position x s and the unperturbed image position Xi, in terms of the Einstein 
angular radius of the isolated BH, 
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where M. = (3.0 ± 0.5) x 10 6 M (Genzel et al. p000|) is the BH mass, Rq = 8.0 ± 0.5 kpc (Reid 
1993 ) is the distance to the GC, D s is the distance to the source, and 6*00 = 1.75" ± 0.20" is the 



Einstein angle for a source at infinity (1" ~ 0.04 pc at the GC). We express the displacement 
vectors ^ and £ p in terms of the Einstein radius of the star, ^/sOe, where e = m*/M, is the mass 
ratio between the star and the BH. This quantity provides a good measure for the scale of the 
stellar lensing zone. Correspondingly, we express areas in the lens plane, c*, in units of eitO'g, the 
Einstein ring area of an isolated star; and express surface number densities of lensing stars in the 
lens plane, £*, by (eftO^) . Areas in the source plane, a s , are expressed in terms of tt6 e , and 
surface number density of sources in the source plane , T, s , are correspondingly expressed in terms 

of HI)- 1 . 

In the limit e <C 1, & is given by the two or four solutions to the equations (Gould &; Loeb 
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The magnification of each of the images is 
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For the GC, we consider solar mass lenses and adopt, for simplicity, a single value of e = 3 x 10 

The angular area in the source plane where a background source would be magnified above a 
given threshold A by the star-BH system, a s (>A), is obtained by identifying the corresponding 
angular area in the lens plane around the unperturbed image and transforming it back to the 
source planeQ . The lensing stars are distributed randomly around the BH and they scan the lens 



1 Source areas that contribute two images which are above the magnification threshold, are counted twice. 



plane as they move. Averaging over time and assuming an isotropic projected stellar distribution 
around the BH, 

1 f dx f 

a s (>A) = - / P(> A,x i9 ) — — dXi. = 2 / P(> A,Xi.)A- l (xi.)xi.dxi. , (5) 

ir J ax im J 

where P(> A,Xi») is the probability for magnification above A, that is, the fraction of dXi, where 
the image is magnified above A. Equivalently, P{> A, X{ 9 ) is the fraction of time a stationary 
source is magnified above A by a closely-passing star. The magnification A, is that due to the BH 

alone, 
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where the positions of these images are related to the source position, 
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Xi.,± = ^\ Xs± V X s + 4 J • ( 7 ) 

De-magnification, i.e. A < A,, is also possible. The probability P is expressed in terms of the 
optical depth 

r* (>A,Xi,) = E*(xj.)cr^(>A,Xj.) , (8) 

where ov(> A, Xj.) is the area in the lens plane where the presence of a star will lead to 
magnification above A, and £*(xi # ) is the stellar lens surface density. In the small optical depth 
limit 

P(> A, Xi.) ~ r*(> A, Xi .) + 6 (A. - A) , (9) 

where is the Heaviside step function. In practice, the lens plane area a+(> A, x-i,) is calculated 
by solving equations (|2])-([|) numerically as functions of $, p in a small test area ao(xi,) around a?j., 
where the effect of the stellar lens magnification is appreciable. The effect of the star outside of uq 
is neglected. With this assumption, the probability for not being affected by any star is exp(— to), 
where To = X*cro- The magnification probability is then given by the sum of the probabilities of 
magnification by a star or by the BH alone, P{> A) = [1 — exp(— T*)] +exp(— tq)Q(A 9 — A), which 
in the limit t*, to <C 1 reduces to equation (J9|), independently of the exact choice of gq. 

For an isolated BH (£* = 0), the corresponding angular area in the high magnification limit, 
A»l, isa s (>A), ^1/2A 2 . 

Equations (|2])-(|]) are valid as long as r^(> A, x^,) <^. 1, i.e. when A is large enough 
so that there is no overlap between the regions of influence (lensing zones) of different stars, 
&+(> A, Xi,) -C S^ 1 (xj.). In addition, these equations ignore the cumulative effect of all the stars 
on the central caustic around the BH (Griest & Safizadeh 1998). For e = 3 x 1CP 7 , the central 
caustic has a negligible contribution to the integrand of equation (|5|) as long as x-i, is not very 
close to unity (i.e. for sources which are not located almost behind the BH). The inclusion of the 
central caustic can only increase our estimated lensing probabilities. 



3. Model 

In §01 below we demonstrate the effect of enhanced microlensing by stars near a massive 
BH by calculating the mean number of lensed images with magnification above a threshold A, 
{Ni(> A)} = cr s (> A)X S , of distant background sources behind SgrA*. To do this, we need to 
specify the projected density of the stellar lenses £*, which determines the optical depth for 
microlensing and hence the cross-section in the source plane a s {> A) (eq. 0), and we also need to 
specify the projected source density T, s and source luminosity function. We express these surface 
densities scaled to #00 = 1.75". 



3.1. Stellar Lenses Around SgrA* 
The stellar mass density distribution near the BH is modeled as a power-law, 

p=(Z-a)p h (pj a , (10) 

with pb = 10 6 Mq pc -3 , r^ = 0.4 pc and an index a in the range of 3/2 to 7/4. The normalization 
of the mass density is based on the dynamic mass measurements of Eckart & Genzel ( |1997[ ) . The 
power-law distribution is indicated by an analysis of the star counts in the inner GC (Alexander 
|1999[ ) and agrees with the theoretical prediction for a relaxed stellar system around a black hole 



(Bahcall & Wolf 1977 ), as is thought to be the case in the GC. For simplicity, we assume that the 
stellar lenses all have the same mass m+ = 1 Mq (e ~ 3 x 10~ 7 ). For these low mass stars, we 
adopt the theoretical prediction of a ~3/2. The corresponding surface mass density (in units of 
M, /nOg) and surface number density (in units of l/e-irOg) at an angular separation x = 6/9e 
from the GC are both given by 

E* = £*x 1 ~ a , (11) 

with X* ~ 0.04 {9e/ &oc) ~ a - The normalization X* is almost independent of a in the above range 
and scales as M" 1 . The ratio of stellar mass enclosed within the Einstein radius and the mass 
of the BH is 2X* /(3 — a) . Note that X* is also the optical depth for lensing by the stars alone, 
neglecting the effect of the BH. As we show below, the actual optical depth increases substantially 
when the shear of the BH is included, an effect analogous to the cross-section enhancement for 
planetary systems (Gould & Loeb 1992). 



3.2. Distant Background Stars 

A rough estimate of the projected density of distant stars, which is based on a model of the 
-fC-band luminosity distribution in the Galaxy (Kent |1992j ) and the assumption that the mean 
stellar AT-band luminosity is solar, suggests that there are of order ~ 100 distant background 
stars within ~ 2" of the BH. We define "distant background stars" as those stars with 6e ^ 1", 



i.e. farther than ~ 4 kpc behind the BH. The surface density of nearby background stars in the 
inner bulge and in the high density central cluster may be as high as ~ 10 3 arcsec -2 (Alexander &: 
Sternberg 1999). However, the Einstein angular radius for stars so close behind the BH is less than 



0.1" and there are not enough lensing stars near the BH on that scale for microlensing to become 
important. Here, we study the role of these close stars as microlenses rather than as sources. 

The stellar K luminosity function (KLF) of the distant background stars can be determined 
from observations in the 6 = 0°,/ = 30° direction (i.e. impact parameter of 4 kpc relative to the 
GC). The infrared Galactic model of Ortiz & Lepine ( |1993j ), which reproduces such observations, 
predicts S s (< 13. 5 m ) ~0.1 (see their Fig. 17b) with a steep slope of b = dlogT, s /dK ~0.4 to 
0.5. The projected stellar population is integrated over most of the Galactic disk and is composed 
of stars of different ages, and so it is reasonable to approximate it as an old, continuously star 



forming population. Observations of such populations in the Galactic Bulge (Tiede et al. |1995| ; 



Holtzman et al. 1998| ) and the central cluster in the GC (Blum, Sellgren & DePoy [1996 ; Davidge 



et al. 1997J ), as well as stellar population synthesis of the GC (Alexander & Sternberg 1999J ) 



indicate that the post-main-sequence KLF, where the ET-luminous sources lie, typically takes the 
form of a single power-law with b ~ 0.3-0.4 down to the main sequence turn-off point at around 
Mkq = 3.3 m , where the power-law flattens and then turns over. Assuming a typical distance of 



2Rq to the distant background sources and Ak ~ 3 m to the GC (Blum et al. |1996[) , and using 
the population synthesis models to estimate the mean stellar K luminosity and the fraction of 
stars less massive than solar, we approximate the KLF as a single power law that extends down 
to a cutoff magnitude of K c ~28.2 mag and normalize it by Kent's integrated luminosity density 
model. The uncertainty in the exact shape of the luminosity function, the patchiness of the 
extinction in the Galactic Plane, and the distribution of star-forming regions in the spiral arms 
behind the GC introduce a large uncertainty to S s . Putting these uncertainties aside, we adopt a 
working model for the KLF of the distant background stars, 

S s (< K) = %1& K , (12) 

with S s = 5 x 10~ 10 (9 e /0og) 2 and b = 0.4. 

The uncertainties inherent in the estimate of T, s are demonstrated by the attempt to 
extrapolate to low luminosity the observed KLF and correct it for distance, projection, foreground 
stars and the higher extinction at / = 0°. The Galactic extinction model of Wainscoat et al. 



( 1992 ) predicts that the integrated K-band extinction in the direction of the GC is ~ 2 m larger 
than that at / = 30°. The correction for foreground stars and for the projection from I = 30° 
to Z = 0° (assuming an exponential Galactic disk with a length scale of 3.5 kpc, Wainscoat et 
al. 1992[) further reduces the projected density by a factor of ~0.5, and the shift from a typical 
distance of ~i? to 2R by another factor of ~ 0.25. This gives S s (< 15.5 m ) ~0.01 (^/floo) 2 for 
the distant background stars, which corresponds to S s = 5 x 10~ 9 (0e/ 9oo) > 10 times higher than 
our adopted normalization (eq. |i~2| ). 

The differential surface density of lensed stars, dT, s /dF\~, is related to their unlensed 



luminosity function by dT, s /dF\~ = A 2 dT ls /dF\~,., where the tilde denotes lensed quantities 
and where dT, s /dF\p oc p- 2 -^- 1 oc jp~ 2 f or the KLF of equation (p^) . Hence, the adopted 
surface density has the critical power-law slope for which the differential flux distribution remains 
unchanged. 



4. Results 

Figure ^a shows the integrated source plane area a s (> A) for producing a lensed image 
of a point source that is magnified by more than A anywhere in the inner 29 e, for 9e = #00 
and the stellar lens surface density model of equation ( |Tl| ) with a = 3/2. The results can be 
scaled to any value of 9e or any stellar lens surface density normalization X* by noting that 
a s {> A) - cr s (> A), oc E* oc 9%~ a (Eqs. §, and @). In particular, for D s = 2R the 



difference in source plane areas is decreased by a factor of 0.6. The highest magnifications are due 
to stellar lenses close to Be- The integration range in equation (|5|) does not include the annulus 
|1 — Xj,| < 5x = 0.005. At that point the semi-analytic approximation (eq. ||]) and the small 
optical depth assumption no longer apply because the change in A, over the angular scale defined 
by a*(> A, Xi,) is no longer small, and r*(>^4, Xj,) <C 1 no longer holds even for very high values 
of A. We verified the validity of our semi-analytic calculations for the case of a BH and a single 
star in the range 5x > 0.001 by comparing cx*(> A, Xi,) as function of A and Xi, with results from 
the direct ray shooting method (Wambsganss & Kubas |2000 ). The two methods give identical 



results to within ~5% down to 5x = 0.1 for both the normalization and slope and remain in good 
qualitative agreement even down to 5x = 0.001 throughout the range where the ray-shooting 
method has sufficiently high resolution. 

For a star at a particular location, the magnification probability approaches the asymptotic 
A~ 2 behavior in accordance with the theoretical asymptotic limit for a point source (eq. [11.21b] 
of Schneider, Ehlers & Falco 1992| ) at some minimum magnification A[,, which depends on the 



external shear at that location. As the projected location of the image approaches 9e, the A~ 2 
tail is recovered at an increasingly larger value of A^. With the above choice for 5x, we find that 
the area a s (>A) initially falls off as A~ 5 ' 2 up to a break at A^ ^ 50 where it begins to converge 
to the asymptotic A~ 2 behavior. In our semi-analytic approach, the value of A\, depends on the 
choice of the truncation radius 5x, with smaller values of 5x corresponding to larger values of A},. 
An exact ray-shooting calculation is required to determine Ab more precisely for a point source. 
Nevertheless, we note that our calculated value for a s (> A) in the asymptotic limit agrees very 
well with the analytic expression for a binary lens (eq. [11.23a] of Schneider, Ehlers & Falco 1992). 



taking into account the fact that our definition of 7 coincides with the shear parameter in the 



limit e -> (Dominik |1999j , Eqs. [12] and [66]) 
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where the integral was evaluated for a = 3/2 on the interval (0 < Xj, < 1 — Sx) U (1 + 6x < Xj, < 2) 
and where E is the complete elliptic integral of the 2nd kind. 

At high magnification the extended nature of the stellar sources can no longer be neglected. 
The maximal magnification for an extended source of angular size 6 S is of order A max ~ W € 1 / 2 8e/0 



(Schneider, Ehlers & Falco 1992 , p. 215), which for a 1 R & star located at 2Rq gives A max ~ 60 
and for a 10 R@ giant at the same location gives A max ~ 20. The absence of the high-magnification 
events is compensated by an increase in a s below A max , which is not taken into account here. We 
find that even for magnifications A <C A mSuX , the cross-section of the BH and stars system is larger 
by factors of ~ 2 for sources at D s = 2Rq, depending on the normalization of the surface density 
model for the lensing stars. 

Figure |2|b shows the differential area da s (> A) /dxj, , for the BH with and without the stars. 
The contribution of stars to the high magnification events extends well beyond the narrow region 
around the Einstein radius that is expected for an isolated BH. 

Figure [||a shows the number of lensed images that will be observed, on average, in the inner 
28 e for different limiting i^-magnitudes, for D s = 2Rq and a = 3/2 
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where we assume that the KLF of equation (|l^) has a sharp cutoff beyond K c , and where 
K A = 2.5 log A, A c = io - 4 ^-^), and a s is parameterized as a s {> A) = d s A~ l - b with a s w 0.16 
(Fig. |^a). An analogous expression with A~ 2 and a lower normalization describes the magnification 
by the BH alone. For low values of A and Kq, (Ni(> A; Kq)) decreases slowly like ai 2, 
because the decrease of a s (> A) is partially offset by the ever larger number of faint stars in the 
population that become accessible. However, once Kq + K A > K c , all the background stellar 
population can be observed when magnified by A, and (iVj(>^4; Kq)) falls off more rapidly, like 
a s (>A)cxA- 1 - 5 . 

The motions of the stellar lenses and the sources introduce time-dependent Poissonian 
fluctuations in iVj(> A; Kq). The fraction of time that at least one lensed image will be observed, 
which is shown in Fig. |b, is P [Ni(> A; K ) > 1] = 1 - exp [- (iVi(>,4; K ))]. 

Present-day observations reach a depth of Kq ~ 16 m (Ghez et al. 1998|) . For this threshold, 
the probability of seeing at least one distant background star in the inner 26 e ~ 2" (D s = 2Rq) 
with A > 5 is small, ~ 0.02%. The lensing statistics will improve significantly with deeper 
observations. We predict that in any snapshot of the central 2" region around the GC to a 
iT-magnitude limit of 23 m , at least one distant background star will be magnified by a factor > 5 
for 10% of the time. 

Variability offers a direct way of identifying the lensed images among the many foreground 



stars. Any source detected with a magnitude Kq at a given position Xj, will be magnified by 
an additional factor of > A for a fraction ~ I /A 2 of the time (Schneider et al. 1992; ). Hence, 



continuous monitoring of detected sources would allow identification of microlensing events. Such 
events can be distinguished from variability of the source stars or from variable patchiness in 
the dust obscuration, through their achromaticity and their generic set of possible time profiles. 
The determination of the microlensing probability and duration distribution can be used to 
constrain the density and velocity distribution of low-mass stars in the GC. We demonstrate 
this connection qualitatively by defining a typical timescale for magnification by more than A as 
t(> A) = ^/ov(> A) I v± oc ■s/m*/ v± (Although the actual event duration depends on the complex 
caustics geometry, we use this definition as a crude estimate). The projected rms transverse 
velocity of the low mass stellar lenses in the GC is (Alexander |1999| ) 

v l ~ 0.37^ , (15) 

V 

where p = 9Rq is the projected distance from the BH. Figure compares t(> A) of magnified 
images of stationary sources lensed by solar mass stars in the shear field of SgrA* with the 
timescale of lensing by isolated stars (with the same velocity field). The duration of a microlensing 
event strongly depend on the position of the image, and is increased from several weeks to several 
months near 9e- The formal divergence of t(> A) near 9e is truncated in practice by the rapidly 
decreasing probability of having a source with an image close to 9e (cf Fig. (2(b) and by the motion 
of the source, which is neglected here. Nevertheless, a sharp increase of the microlensing timescales 
is expected near 9e- In contrast with timescales of lensing by isolated stars, the timescales inside 
9e fall rapidly to zero in the presence of the BH because the strong de-magnification by the BH 
(eq. 0) has to be compensated by very high stellar lens magnification. The lensed lightcurves 
should show complex structures, unlike the symmetric smooth lightcurves of an isolated point 



mass lens. Examples of such lightcurves can be found in the literature (Mao & Paczynski 1991 



Gould & Loeb 1992; Bolatto & Falco 1994, and in particular Wambsganss 1997). 



5. Discussion and Summary 

The presence of a dense stellar cluster around a super-massive BH changes considerably the 
lensing properties of both the stars and the BH: the magnification cross-section, the magnification 
probability, the variability timescale and the lightcurve structure. The modification of the 
microlensing properties of the stars is most pronounced near the Einstein radius of the BH. 
Previous work on lensing by the BH alone (Wardle & Yusef-Zadeh 1992| ; Alexander & Sternberg 



1999 ) focused on stellar sources from the inner central cluster, for which 9e ~ 0.02". Here we have 
considered lensing of distant background stars by the BH and the stars around it, for which the 
angular scale is 9e ~ 2". 

We applied our calculations to the BH and the dense central cluster in the GC, and found 
that the lensing probability is increased over that of an isolated BH by factors ranging from ~ 2 



10 



for low magnifications (;> 5) to ~ 3 for high magnifications (~ 50). We estimate that in any 
snapshot of the central 2" region around the GC to a .ff -magnitude limit of 23 m , at least one 
distant background star will be magnified by a factor ;> 5 for 10% of the time. The duration of a 
microlensing event can be lengthened by an order of magnitude near the Einstein radius of the BH. 
Moreover, the gravitational shear of the BH changes considerably the shape of the microlensing 
lightcurves so that they are no longer symmetric and smooth as they are for isolated stars. 

Variability provides the distinguishing signature of microlensed stars. Image subtraction can 
be used to pick out lensed images of distant background sources from the many foreground stars. 
The achromaticity of the lightcurves and their generic (non-periodic) temporal structures can be 
used to separate them from those of intrinsically variable stars. Since the variability timescale 
scales as yjm+/v_\_, the lightcurves contain information on the stellar mass function, the stellar 
density, and the stellar velocity dispersion as functions of the projected distance from the BH. In 
particular, the lensing timescales could probe mass segregation deep in the potential of the BH 



(Miralda-Escude & Gould 2000). The vast majority of the stellar lenses are low-mass faint stars 
that are many magnitudes dimmer than the current detection limit. Microlensing can probe the 
statistical properties of these yet unobserved stars. 

It is important to note that our quantitative results are subject to a number of approximations 
and uncertainties. For the sake of simplicity, we ignored the fact that the background stars may 
have a broad distribution of distances and hence different values for the Einstein radius. In 
addition, the model for the surface density and luminosity function of the distant background stars 
suffers from large uncertainties given the quality of present-day data. We did not attempt to deal 
in detail with the finite size of the stellar sources, which becomes relevant at high magnifications. 
We also ignored the contribution of blended light from the lensing star to the inferred flux from 
the source. However, this contribution is likely to be small for bright sources and low-mass lenses. 
Lastly, the semi-analytical formalism for calculating the magnification is not valid very near 9e 
and so our quantitative predictions for {Ni{> A)) are less reliable for very high values of A. More 



precise calculations, for example by the ray-shooting method (Wambsganss 1997; Wambsganss & 



Kubas 2000 ) will be required to improve the accuracy of the predictions. 



Finally, we emphasize that our calculations can be applied to any extragalactic supermassive 
black hole. For example, the core of M87 harbors a black hole mass of ~ 2 x 1O 9 M . The Einstein 
radius of this black hole for sources at infinity is ~ 1", close to that of SgrA*, and there is no sign 
for obscuration by dust within its boundary. However, it is more likely that the sources will be 
stars in M87 a few kpc behind the black hole, for which 6^ ~ 0.03". For the higher black hole 
mass in M87, e ~ 5 x 10~ 10 and so microlensing events of stellar sources which are embedded in 
that galaxy will have durations comparable to those around SgrA* since t ~ (e/M,) ' (Di6 e ) ' , 
where D\ is the distance from the observer to the BH. Since crowding and flux sensitivity severely 
limit the detection of individual stars at the distance of M87, one may search for microlensing 
of surface brightness fluctuations, the so-called pixel lensing (see, e.g. Gould 1996). The Next 
Generation Space Telescope (NGST; [http://ngst.gsfc.nasa. gov/1) , with its sub-nJy sensitivity in 
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the wavelength band of ~ 1-3. 5^m and its ~ 0.06" resolution, is ideally suited for such a study. 

We thank Joachim Wambsganss for communicating results from ray-shooting calculations 
prior to publication (Wambsganss & Kubas 2000) and an anonymous referee for useful comments. 
The KLF source model was revised following results communicated to us by J. Chaname, A. 
Gould & J. Miralda-Escude prior to publication. AL thanks the Institute for Advanced Study for 
its kind hospitality during the initiation of this paper. This work was supported in part by NASA 
grants NAG 5-7768, 5-7039 and NSF grants AST-0071019, AST-0071019 (for AL). 
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Source plane 



1st unperturbed image 




2nd unperturbed image 



Perturbed image 



Observer 



Fig. 1. — A sketch defining the notation used in this paper. The presence of a perturbing star at 
position £ p relative to the unperturbed image at Xi, splits this image into 2 or 4 images at positions 
£j relative to the unperturbed image. For clarity, only one of the multiple images due the the star 
is shown and the proportions are exaggerated. 
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Fig. 2. — (a) The mean area in the source plane for magnification by more than A by the BH 
alone, and by the BH and stars of an image in the inner 29 e, for 9e = #oo and the stellar lens 
surface density model of equation (|ll|) with a = 3/2 (see text for scaling to arbitrary 9e)- The 
two vertical hash marks roughly indicate the maximal magnification possible for extended sources, 
a 1 Rq and a 10 Rq star at a distance of 2Rq. (b) The differential source area per annulus in the 
lens plane that is magnified by more than A. The case of a BH and stars (full lines) is compared 
to the contribution from the BH alone (dashed lines) for three threshold values of A. 
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Fig. 3. — (a) The average number of lensed images magnified by more than A that will be observed 
in the inner 29 e with a limiting K-band magnitude Kq, for D s = 2Rq and the stellar lens surface 
density model of equation (11) with a = 3/2. (b) The fraction of time that at least one lensed image 
magnified by more than A, will be observed in the inner 28e with a limiting K-band magnitude Kq 
(Pel scales linearly with T, s ). The labeled solid lines are for the BH and stars and the unlabeled 
dotted lines, which follow the same Kq order, are for the BH alone. 



16 



300 



250 



200 



>, 

CO 
"D 



150 



100 



50 



J 


A>5 
A>10 
A>20 I 






-^^=^Ssi^T^ 


1 — " 



1.5 



2.5 



6 (arcsec) 



Fig. 4. — The typical duration of a microlensing event for a stationary source at infinity by a solar 
mass lens (assuming the mean rms velocity of equation ]l5|), as a function of angular separation 
9 = Xi,6 'oo from SgrA*. Curves are shown for different magnification values A, comparing stars in 
the shear field of the BH (full lines) with isolated stars moving at the same velocity (dashed lines) . 



